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Abstract 

The Rancho La Brea Tar Pits is the world's richest and most important Late Pleistocene fossil locality and best renowned for 
numerous fossil mammals and birds excavated over the past century. Less researched are insects, even though these 
specimens frequently serve as the most valuable paleoenvironemental indicators due to their narrow climate restrictions 
and life cycles. Our goal was to examine fossil material that included insect-plant associations, and thus an even higher 
potential for significant paleoenviromental data. Micro-CT scans of two exceptionally preserved leafcutter bee nest cells 
from the Rancho La Brea Tar Pits in Los Angeles, California reveal intact pupae dated between ~23,000-40,000 radiocarbon 
years BP. Here identified as best matched to Megachile {Litomegachile) gentilis Cresson (Hymenoptera: Megachilidae) based 
on environmental niche models as well as morphometries, the nest cells (LACMRLP 388E) document rare preservation and 
life-stage. The result of complex plant-insect interactions, they offer new insights into the environment of the Late 
Pleistocene in southern California. The remarkable preservation of the nest cells suggests they were assembled and nested 
in the ground where they were excavated. The four different types of dicotyledonous leaves used to construct the cells were 
likely collected in close proximity to the nest and infer a wooded or riparian habitat with sufficient pollen sources for larval 
provisions. LACMRLP 388E is the first record of fossil Megachile Latreille cells with pupae. Consequently, it provides a pre- 
modern age location for a Nearctic group, whose phylogenetic relationships and biogeographic history remain poorly 
understood. Megachile gentilis appears to respond to climate change as it has expanded its distribution across elevation 
gradients over time as estimated by habitat suitability comparisons between low and high elevations; it currently inhabits 
mesic habitats which occurred at a lower elevation during the Last Glacial Maximum ~21,000 years ago. Nevertheless, the 
broad ecological niche of M. gentilis appears to have remained stable. 
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introduction 

Megachile Latreille [1] is a large, worldwide genus of approxi- 
mately 1,500 species of largely leafcutting, solitary bees. In the 
Western Hemisphere they inhabit temperate, arid, and tropical 
regions extending from Alaska to Tierra del Fuego [2] . There are 
1 1 8 species native to North America [2] . The abundance of 
megachilids in California is not surprising given the wide diversity 
of habitats and microclimates [3,4]. 

Leafcutter bees are named for their use of leaf pieces in nest 
building. They constitute Megachile belonging to Michener's Group 
1 [5] in which bees frequently construct two or more cells in a 
linear series. Their nesting sites are found under the bark of dead 
trees, in stems, in the burrows of wood-boring insects or in 
burrows self-dug in loose soil or those made by other animals 
[2,5,6]. The females use their sharp, serrated, scissor-like 



mandibles to cut oblong and circular leaf pieces, most likely from 
plant sources near the nest [2,5]. They line the nest cavities with 
overlapping layers of the oblong-shaped leaf disks. The leaf edges 
are compressed to extrude sap that, in combination with saliva, 
creates a glue-like substance that keeps the cells sealed and intact 
[2] . Each cell is provisioned with pollen and nectar by the female 
before she deposits a single egg on the food mass. After depositing 
the egg, she seals the cell with one to several circular leaf "caps" 
[2,5,6] . After a few weeks, depending upon species, the eggs hatch, 
and the larvae develop through multiple instars and feed on the 
provisions. Mature larvae spin cocoons of two or more layers of 
silk and diapause as prepupae. Cocoons are sturdy structures [7] 
made increasingly airtight by the larva's secretion of a brown 
liquid that fdls and hardens the interstices of the silk layers [6]. 
This application binds the silk mesh and makes the cocoon 
extremely durable. Simultaneously, this fastens the cocoon's outer 
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surface to the surrounding leaf disks that firmly hold the structure 
together. The larvae subsequently pupate and emerge as adults by 
chewing their way out through the cap. 

That females may spend the majority of their time collecting 
pollen and nectar to provision their young [2] and construct 
intricate nests with specific materials indicates a very complex and 
highly evolved plant-insect interaction, and strongly suggests a 
long evolutionary history [8]. The use of leaf disks of various sizes, 
shapes, and textures also reflects highly complex and evolved 
behavior [2,8,9,10]. 

As currendy known, the megachUid fossil record is restricted to 
the Cenozoic based on body fossils preserved as compressions and 
three-dimensionaUy preserved in amber, as well as trace evidence 
from fossil angiosperm leaves whose margins show smooth-edged 
oblong and circular cutouts [8,11,12]. Engel [13-15] and Engel 
and Perkovsky [1 1] have compiled the evolutionary history and an 
over\dew of the body fossil record, respectively. Morphological 
data (body fossils and leaf cutouts attributed to Megachile) and 
molecular data do not always agree on the time divergence of the 
genus [8,11-14,16-22]. Although the phylogenetic relationships 
and evolutionary history of the genus have become clarified as 
more studies incorporate molecular data [23], the fossil record 
remains incomplete and some specimens assigned to MegachUidae 
may need revision. For example, molecular data suggest that 
Megachihdae arose in the Cretaceous about 140-100 mya, but the 
genus Megachik is estimated to have originated only 22 mya [23] . 
However, leafcutters are derived species of Megachik and therefore, 
the fossil record based on leaf cutouts from the Early to middle 
Eocene in North American and Europe [16-20] suggests that 
basal divergences in the MegachHini occurred earlier in the 
Paleocene or Latest Cretaceous [11]. 

Here we report on fossil Megachile nest ceUs with pupae 
(LACMRLP 388E) recovered from the Late Pleistocene Rancho 
La Brea Tar Pits in southern California. Though geologically 
young, this is the first report of three-dimensionaUy preserved 
Megachik nest cells that shows rare preservation and life-stages. 
The pupal morphology, nesting behavior, and cell construction of 
LACMRLP 388E best match M. gentilis, a member of the native 
Nearctic Litomegachik MitcheU [24], a subgenus that today ranges 
from southern Canada to Cuba and southern Mexico [5,25,26]. 

The discovery of LACMRLP 388E provides valuable informa- 
tion for better understanding the environmental conditions of 
southern California in the Late Pleistocene. By setting specimens 
within a geological as wefl as an ecological context, Quaternary 
fossils are shown to be valuable precursors to modern biota [27]. 
Although the asphaltic deposits at Rancho La Brea are most often 
associated with vertebrate remains from saber-toothed cats and 
mammoths, the insects and plants found there are also significant 
fossils because they are original material, and thus, intact, three- 
dimensional, and structurally complex. As such they can provide 
the most valuable paleoenvironmental information for the richest 
Ice Age fossil locality. Our goal was to syn[:hronize data by 
identifying both nest cell insect and plant material in order to 
make the significant paleoenvironmental inferences possible. 
These elforts resulted in environmental data, that a mesic 
environment occurred at a lower elevation than today at Rancho 
La Brea, the well-estabUshed provenance for LACMRLP 388E. 
This research also resulted in new information on M. gentilis, 
including diagnostic features of its nest ceU architecture and insight 
regarding its relatively conserved ecological niche since the Late 
Pleistocene. 

In additional to its role as a sensitive paleoenvironemental 
indicator and providing new information on M. gentilis, 
LACMRLP 388E is of rare value because its remarkable 



preservation, especially of the intact pupae, is of a standard 
unusual even for fossil material from Rancho La Brea. 

Results 

LACMRLP 388E 

Provenance and age. LACMRLP 388E was excavated in 
1970 from the northeast corner of Pit 91, one of approximately 
100 asphaltic deposits at Rancho La Brea in Los Angeles. It was 
found as one intact specimen (LACMRLP 388E) in an area with 
asphalt-impregnated sediment and bone at a depth of 205 cm. 
However, on subsequent handling the two nest ceUs separated 
(LACMRLP 388Ea, Fig. \A-C; LACMRLP 388Eb, Fig. \D-G,J) 
along with a portion of the outer leaf layers (LACMRLP 388Ec) 
(Fig. \H, 7). 

Fifty fossil bones from Pit 91 have been dated between 23,000- 
40,000 years BP based on radiocarbon dated bone collagen [28]. 
This is with the exception of one, anomalous bone on the opposite 
side of the pit from LACMRLP 388E, which was dated to 14,000 
years BP [28]. Furthermore, LACMRLP 388E was recovered 
from a depth of less than 243 centimeters that included a cluster of 
bones dated to ages older than 35,000 years BP, increasing the 
likelihood that LACMRLP 388E is between 35,000 and 40,000 
radiocarbon years BP [28]. That the nest cells were constructed 
underground but close to the surface before becoming embalmed 
in asphalt-rich matrix indicates that they may have been near an 
asphalt pipe from which oil soaked into the surrounding sediment 
[28]. 

Many aspects of LACMRLP 388E nest cell construction 
increased the probability of fossUization in asphalt. An under- 
ground burrow provided a protected environment from the 
elements, as well as many predators. Glue-like saliva containing 
anti-bacterial and anti-fungal properties likely prevented decay 
before fossflization [5,29,30]. Water-resistance was created from 
sap extruded from the ceU leaf pieces, tight assembly of multiple 
leaf layers, the hydroscopic properties of the leaf epidermis, and 
the internal silk cocoon which was hardened and airtight from 
larval secretions plus lipids from the pollen provisions. Finally, the 
life stage-a bee in the pupal stage-would be less vulnerable to 
desiccation than a larva whose exoskeleton is less strongly 
sclerotized. 

Pupae. Originally identified on-site as a bud (Obermayr, Pit 
91 field notes p. 1770, 1970), LACMRLP 388E was later 
tentatively identified to the bee family Apidae. Micro-CT scans 
of LACMRLP 388E reveal a pupa within each cell, indicating that 
the specimen was preserved sometime during flowering season 
(Fig. 2A, C and Fig. 3A, C, Videos SI -6) since Megachik do not 
overwinter as pupae. A comparison of LACMRLP 388E to a 
modern leafcutter bee {Megachile rotundata) show obvious similarity 
in general morphology at pupal stage (Fig. 4). While nest ceUs are 
designed to house Megachik offspring, the megachHid cleptopar- 
asite Coelioxys may usurp nest ceUs. Images of the pupae in dorsal 
and ventral views (Fig. 2A, B, F and Fig. M, B, F, Videos S3, 4) 
have the same general body morphology, presence of setal pits on 
scutum, the long labrum, and oval shape of the metasoma of 
Megachik, confirming the identity of the nests based on ceU 
construction and structure. The shape of the metasoma and 
posteriorly rounded mesosoma eliminates the possibility of 
Coelioxys, which would have a more elongate, tapered metasoma 
in dorsal and ventral views and pointed axillae. One pupa is 
present in each ceU and both sexes are represented (Fig. 2 and 
Fig. 3, Videos Sl-6), the male being distinguished by a dorsal 
horizontal ridge on metasomal segment six (Fig. 2. A, B). Because 
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Figure 1 . Images showing nest cell construction and modern nest cell of Megachile gentilisiot comparison. A-C = LACMRLP 388Ea; D-G, 
J = LACMRLP 388Eb; H, I = LACMRLP 388Ec. (A) Nest cell containing male pupa showing cylindrical shape, tapered, rounded bottom at left typical of 
the first constructed cell, and remains of oblong leaf disc with Type 1 venation. (B) Bottom of first constructed cell (containing male pupa) with 
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possible portion of bottom circular leaf disk visible and outlined with arrows. (C) Cap of nest cell containing male pupa. (D) Nest cell containing 
female pupa. Arrow shows margin of oblong leaf disc with Type 2 venation. (E) Circular, bottom disc of nest cell (containing female pupa). In life, this 
end abutted the anterior end of LACMRLP 388Ea. (F) Cap of nest cell containing female pupa. (G) Nest cell of female pupa showing oblong side wall 
leaf cutout which does not reach bottom of cell and is instead supported by circular bottom disc. (H) Remains of oblong leaf disc with relatively 
smooth-cut margins and Type 3 venation. (I) View showing five overlapping oblong disks (1-5) comprising the sidewalls and circular bottom disk, 
discs. (J) Nest cell containing female pupa showing Type 4 venation on upper, right corner (K) Nest cell of modern M. gentiUs, showing circular disc 
bottom and oblong, sidewall leaf 
doi:10.1371/journal.pone.0094724.g001 



the cells were connected in situ, and are not cleptoparasitized, they 
can only be the same species. 

Nest Cell Construction 

Nest cell morphology. The two nest cells that constitute 
LACMRLP 388E are approximately 10.5 mm in length and 
4.9 mm in diameter, constructed of oblong and circular leaf 
cutouts (Fig. lA-J). The nest cells when found were joined end to 
end, the male cell in Figure \A being the cell first assembled in the 



ground, while the female cell (Fig. \D, G, J) would have been 
closer to the nest entrance. They were held together by additional 
cutouts that lined the nest cells but some dislodged after excavation 
and were lost. With subsequent handling another small portion 
sloughed off but was conserved. It consists of an additional five 
oblong-shaped disks (Fig. \H, I). Outward morphology (Fig. \A-G, 
J) and cross-sectional view of the cells (Fig. 2E and Fig. 3-E, Video 
S7) show the sides and inner lining consist of 5-7 overlapping 
smooth-edged, oblong-shaped disks, the ends capped by small 




Figure 2. Micro-CT scans of LACMRLP 38SEa showing male pupa and its position within the nest cell. (A) Dorsal view of pupa within 
nest. (B) Dorsal view of pupa. (C) Lateral view of pupa within nest. (D) Lateral view of pupa. (E) Cross-section of nest and pupae. (F) Ventral view of 
pupa. 

doi:1 0.1 371 /journal.pone.0094724.g002 
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Figure 3. Micro-CT scans of LACMRLP 388Eb showing female pupa and its position witKiin tKie nest cell. (A) Dorsal view of pupa within 
nest. (B) Dorsal view of pupa. (C) Lateral view of pupa within nest. (D) Lateral view of pupa. (E) Cross-section of nest and pupae. (F) Ventral view of 
pupa. 

doi:1 0.1 371 /journal.pone.0094724.g003 



circular disks (Fig. IC, F, Video S7). At least three discs cap the 
end of tlie male cell (Fig. li^. 

Most often, oblong leaves form the side walls and are bent into a 
cup at the proximal end (furthest from the nest entrance) that is 
glued with saliva and extruded leaf sap. A cap consisting of a series 
of layered circular discs is placed on the distal end (closest to the 
nest entrance). However, the separation of the two cells reveals the 
nest cell containing the female pupa includes a circular cap 
(Fig. IE), as well as a circular bottom (Fig. \F). Certain Megachile 
species insert a circular disc to provide an internal brace for the 
nest cell if the cup pieces are looser and less glued together than as 
constructed in other species [9,10,31]. Therefore, we presume that 
the nest cell containing the male pupa was the first cell in the 
cavity because it is hard to distinguish if the circular base is formed 
from the folded edges of the larger leaf fragments or the insertion 
of a smaller internal disc (Fig. I A, B arrows). Literature on 
Megachile nest cell construction informs that the inclusion of a 
circular base is an uncommon construction for the genus [2,5]. 



Among extant Nearctic Megachile the small size of the cells and 
the use of leaves in nest construction limit candidate species to 
small members of the subgenus Litomegachile, specifically M. brevis 
Say, M. gentilis (Fig. \Kj, M. onobtychidis CockereU and M. 
pseudobrevis Mitchell. Megachile brevis and M. pseudobrevis are excluded 
from consideration because the former uses petals in cell 
construction [32] and the latter's current distribution is restricted 
to a limited range in the southeastern United States [25]. 
Exclusion of M. pseudobrevis is further supported by the incorpo- 
ration of petals in two nests reported from Florida [33]. Further 
examination of the distinct nest cell architecture and exclusive use 
of leaves as ceU material for LACMRLP 388E (Video S7), in 
combination with pupal morphology (Videos SI -6), indicate that 
these cells belong to M. gentilis, a species whose current distribution 
is concentrated in the southern portion of the Pacific coast, the 
Southwest, northern Mexico, and Texas [25, Discover Life 
Megachile gentilis, http://www.discoverlife.org/mp/ 
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20m?kind = Megachile+gentilis, 2013] rather than M. onobrychidis 
(see below). 

Plant material used. Since Megachile typically cut segments 
from the leaf margin, we can identify at least four, possibly five 
different taxa used in construction of the two cells based on the 
dilferent patterns of venation (Figs. \A, D, G, J). Dicotyledonous 
leaves typically have one or more large primary veins, with a 
reticulate network of successively narrower second, third, fourth, 
and fifth order veins in a pattern consistent within a species [34] . 
Although leaf venation is especially useful in fossil plant 
identification, venation chara(ters are most reliable in combina- 
tion with other morphological details such as leaf shape and size, 
shape of the leaf apex and base, and margin and petiole features 
[34], characters unavailable in the small areas of the leaf cutouts. 
Regrettably, this makes it difficult to identify with confidence the 
fragments to a specific species or even family. Howcvc-r, leaf 
texture, thickness, and the venation that is present are consistent 
with the cutouts being from woody trees, shrubs, or vines rather 
than herbaceous taxa, the leaves of which do not preserve well as 
fossils. Some Megachik include petals in their cells but venation of 
the cutouts of LACMRLP 388E is not consistent with petal 
venation [35-37], nor is their thickness (Video S7). 

Historic and Contemporary Environmental Niche Space 

We further researched the best species match for the La Brea 
bees by examining the historic environmental niche space of the 
La Brea Litomegachile relative to the contemporary environmental 
niche of M. gentilis and M. onohychidis, the species whose nest cell 
morphology, distribution, size, and nesting behavior is closest to 
M. gentilis. To estimate the contemporary environmental niche 
space of M. gentilis and M. onobrychidis we associated a suite of 
informative bioclimatic variables with georeferenced natural 
history records (see Materials and Methods). We then constructed 
the historic environmental niche space of the La Brea bees by 
associating bioclimatic variables values of the Last Glacial 
Maximum (LGM, ~ 2 1,000 years ago) with georeferenced 
coordinates of the Rancho La Brea Tar Pits (latitude = 34.063068, 
longitude =— 118.355412). While both species are broadly sym- 
patric throughout California at present (Fig. 5), we found that the 
historic environmental niche of the La Brea bees in the Los 
Angeles Basin during the LGM fell within the 95% data 
concentration ellipse of the contemporary environmental niche 
space for M. gentilis in a principal component (PC) analysis (PC 
1=32.7%, PC 2 = 26.5%) (Fig. 6). Furthermore, LGM habitat 
suitability estimates based on contemporary habitat suitability 
models (HSM) suggest that M. gentilis had a higher probability of 
being distributed at low elevations in the Los Angeles Basin than 
M. onobrychidis (Fig. 5, Fig. S2) (see Material and Methods). 

Due to broad overlap in geography and environmental niche 
space between M. gentilis and Af. onobrychidis, we further tested the 
assignment of the historic environmental niche of the La Brea bees 
using cluster analysis. Ninety percent of the contemporary M. 
gentilis distribution records, along with the La Brea fossU 
distribution record and associated historic environmental space 
values, were assigned to the contemporary environmental niche 
space of Af. gentilis (Fig. 6). Sixty-five percent of the contemporary 
M. onobrychidis distribution records were assigned to the contem- 
porary environmental niche space of M. onobrychidis (Fig. 6). The 
high assignment of the contemporary M. gentilis records to the M. 
gentilis environmental niche space is likely due to the narrow 
environmental distribution of the species. In comparison to M. 
onobrychidis, M. gentilis is distributed across mesic environments, 
where temperatures rarely fall below freezing during the coldest 
months (Fig. SI), with ample precipitation during the wettest 



months (Fig. SI). The ability for M. onobrychidis to exist in highly 
variable environments may have allowed that species to endure 
prolonged, freezing temperatures during overwintering periods, 
and disperse across a broad range of environments, whereas the 
environmental specificity of M. gentilis supports its narrow 
distribution in the Pacific US (Fig. 5). Cluster and principal 
component analysis, as well as estimates of habitat suitability 
during the LGM are consistent with our species identification of 
the La Brea bees unearthed in the Rancho La Brea Tar Pits in 
southern California as M. gentilis. 

Discussion 

The provenance of specimens recovered from Rancho La Brea 
cannot always be discerned because alluvial wash may have 
deposited some of the recovered material at a distance from their 
original location [38]. Nevertheless, the remarkable preservation 
and lack of water damage to LACMRLP 388E indicates with near 
certainty that it was assembled in the ground on the site of Pit 9 1 
where it was found. The provenance of the specimen, in 
combination \\ ith emironmental niche models and morphometric 
investigations, constitutes strong and comprehensive evidence to 
support confident fossil identifications and subsequent paleoenvi- 
ronmental inferences for the Late Pleistocene, an epoch of 
particular climatic variation. 

Because the leaves used in nest construction were collected in 
close proximity to each other and at the site of deposition, we can 
infer from their identification that a woody and riparian habitat 
existed at Rancho La Brea ~23,000-40,000 radiocarbon years 
BP. This inference is supported by Maxent habitat suitability 
models which suggest that M. gentilis was distributed in a mesic 
environment that likely occurred at a lower elevation during the 
Last Glai:ial Maximum (LGM) -21,000 radioi:arl)()n )'(;ars BP 
(Fig. SI, Fig. S2). Presently, M. gentilis occurs in higher elevation 
habitats that constitute the perimeter of its projected habitat 
suitability during the LGM (Fig. 4). These areas surrounds the Los 
Angeles Basin in which Rancho La Brea is situated (Fig. 4). 

Fossil plants from Pit 91 excavations have been grouped into 
four categories including a riparian association which includes 
plants from flood plain and canyon habitats [39]. Aquatic and 
moisture-loving plants indicate a permanent water source near 
Rancho La Brea [39,40]. But additional plants from Pit 91 
indicate other associations such as coastal sage scrub, chaparral, 
and deep canyon [39-41]. Most of the excavated plants probably 
inhabited Rancho La Brea, but others may have been transported 
by streams or floods [39,41]. LACMRLP 388E suggests for die 
immediate Rancho La Brea area the presence of a gallery forest 
within a riparian zone. These forests exist along watercourses often 
in arid regions and standout in contrast to e.g. an adjoining 
grassland or open woodland habitat. The indication of a mesic 
environment at Rancho La Brea is relatively consistent with 
coastal run-olf records [42] which indicate increased precipitation, 
possibly associated with increased glaciation [43], at ~25,000- 
20,000 years ago. Relative comparison of the geographic 
distribution of habitat suitability between tlu; contemporary and 
LGM models reveals that M. gentilis has remained virtually unique 
to California and southwestern Arizona (Fig. 5) where the average 
temperature of the coldest quarter rarely falls below freezing (0°C). 
The contemporary distribution of M. gentilis is limited to mesic 
environments where the diflerence between the maximum of the 
warmest month and the minimum of the coldest month is centered 
around 33°C (Temperature Annual Range, Fig. SI). That Ad. 
gentilis has retained its abiotic niche since the LGM provides a 
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Figure 4. Comparison between (A) a female modern leafcutter bee (Megachile rotundata) and (B) LACMRLP 388E male and in order 
to display obvious similarity in general morphology at pupal stage. 

doi:1 0.1 371 /joumal.pone.0094724.g004 



strong benchmark with which to compare the climate restrictions 
of other fossils from Rancho La Brea. 

Most animals and plants excavated from Rancho La Brea are 
extant, and those that have become extinct are mostly mammals 
such as common, larger carnivores (e.g. saber-toothed cats and 
dire wolves), common, medium-to-large herbivores (e.g. the 
western horse, mastodonts and mammoths), and some birds 
[39]. Only two species of insects recovered from Rancho La Brea, 
scarabs that relied on dung from mammals that became extinct, 
may have died out as well [39,44]. The current geographic 
distribution of most of the insect species from Rancho La Brea 
occurs in warmer parts of the United States and Mexico. Insect 
damage to the fossil bones attributed to tenebrionid and dermestid 
beeties is consistent with a warm period of at least 4—5 months at 
the time that the fossils were trapped [45]. This differs from 
interpretations of Late Pleistocene climatic conditions in southern 



California based on poUen from deep sea cores which indicate a 
shift to a cooler environment, as indicated by a transition from 
predominantly hardwood scrub-oak vegetation to coniferous pines 
and juniper [46-48] during the Last Glacial Interval from 24,000— 
14,000 years BP [46]. However, offshore palynological records 
may not reveal short-term increases in temperature that would 
result in both the colonization of insects with warm climate 
restrictions for brief periods due to their higher mobility, and in 
increased asphalt entrapment from more active seeps [45]. The 
current geographic distribution and climate restrictions of insect 
species preserved at Rancho La Brea, along with the physical 
properties of natural asphalt, indicate that the fossil insects 
represent warmer intervals of the Late Pleistocene [45] . But this 
does not resolve the discrepancy that some fossils recovered from 
Rancho La Brea such as M. gentilis inhabit mesic environments, 
while most of the tenebrionid beetles preserved at Rancho La Brea 
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Figure 5. Estimated geographic distribution of M. gentilis anA M. onobrychidis a\ present and during tKie last glacial maximum (LGM; 
—21 000 years before present) in the western United States. Darker colors represent high MaxEnt habitat suitability (values closer to H = 1), 
whereas lighter colors represent low HS (values closer to HS = 0). Black squares represent georeferenced point localities of the respective species 
(Table SI). Red crosshairs represent the location of the La Brea Tar Pits in Los Angeles, California. Habitat suitability distributions of both species are 
based on eight bioclimatic variables. 
doi:1 0.1 371 /journal.pone.0094724.g005 



inhabit dry scrub and woodland habitats [49] . Further research on 
insects that can be dated to a specific time period, and that also 
have a clear provenance for Rancho La Brea, as well as narrow 
climate restrictions, will greatly enrich our understanding of the 
paleoenvironment, as well as climate variation during the Late 
Pleistocene in southern California. 

Aside from providing paleoenvironmental data, LACMRLP 
388E enhances the fossil record of megachUids and the genus 
Megachile by providing the first reported specimen from the 
Pleistocene. This contributes to a better understanding of a genus 
whose phylogenetic relationships are still being researched. 

The examination of LACMRLP 388E also presents neontolo- 
gists with important information. O'Toole and Raw [2] state that, 
despite the diversity of nest sites and habitats exploited by 
leafcutter bees, their nest cells "dilfer litde in structure and nesting 
behavior," a trait that makes the cells often difficult to identify. Yet 
the nest cell morphology of LACMRLP 388E provides one of the 
most important diagnostic tools to species-level. Nest cell 



construction of leafcutter bees is rarely closely examined or 
recorded. Where this has occurred [9,10,31], species-specific 
aspects of nest cell architecture have been documented. We 
surmise that there is indeed great morphological variation among 
Megachile species as has been demonstrated for another megachHid 
genus, Osmia [50] . A record of species-specific materials and nest 
structure may provide diagnostic tools useful in the field and for an 
organism that is often difficult to identify at an immature stage. 

Materials and Methods 

Excavation and Preparation 

The cells, which were nested in chunky asphalt-impregnated 
matrix, were easier to excavate than if they had been preserved in 
more asphalt-concentrated matrix and therefore, were not exposed 
to solvent used in other preparation techniques. In addition, the 
immediate coating of the specimens with glyptol, an alkyd resin, 
upon discovery may have prevented further breakage during 
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Figure 6. Principal components analysis of environmental space use by M. gentilis and M. onobrychidis. Estimates of environmental space 
use include eight bioclimatic variables, as well as, Cartesian {i.e., latitude and longitude) and elevation variables. Clusters reflect 99% data 
concentration ellipses determined with the scatterplot (car library) functions in R v3.3.0 (R Core Team 2012). 
doi:10.1371/journal.pone.0094724.g006 



handling and preparation (Obermayr, Pit 91 field notes p. 1770, 
1970). 

Micro-CT Scanning and Reconstruction 

A 3D model was built for each pupae using Micro-CT data 
(Videos Sl-6). Specimens were scanned on a SCANCO Micro-CT 
Scanner Model VI. 2a. Linear Attenuation was 1/cm, 7.000000e+ 
01 kVp with a slice thickness of 5 microns. Micro-CT slices were 
analyzed and 3D reconstructions of each bee buUt using Amira 
5.4. LACMRLP 388Ea (male) was output as 2172 slices and 
reconstructed using 3,021,012 triangles (Videos S3, 5). 



LACMRLP 388Eb (female) was output as 2849 sUces and 
reconstructed using 4,380,432 triangles (Videos S4, 6). To build 
3D models of the nests, we scanned each nest on a MicroCAT II 
small animal CT system (Siemens Preclinical Solutions, Knoxville, 
TN, USA). Exposure settings were 70 kVp, 500 niAs, with a slice 
thickness of 2 mm and output as 436 slices. The female nest was 
buUt and output at 18,540 triangles. The male nest was built and 
output at 1 1,012 triangles. AH models were cleaned and smoothed 
in Geomagic 10. 
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Use of Comparative Material 

The fossil specimens were also compared to material from the 
comprehensive collection at the USDA-ARS Pollinating Insect 
Research Lab. The plant materials used and specific architecture 
was examined from the nest cells of other Megachile species of 
similar, diminutive size that include western United States 
distribution. 

Estimation of Habitat Suitability at Present and During 
the LGM 

Environmental space use. To investigate differences in 
environmental space use between M. gentilis and M. onobrychidis, we 
performed principal components analysis of bioclimatic variables, 
altitude, and Cartesian coordinates [i.e., decimal degrees latitude 
and longitude) associated with distribution records. Contemporary 
environmental space use was estimated from recent specimen 
locality records (Dataset SI; National Pollinating Insect Database, 
Logan, UT) for eight informative bioclimatic variables: mean 
temperature of warmest quarter, mean temperature of coldest 
quarter, mean temperature of wettest quarter, mean temperature 
of driest quarter, precipitation of \\Tttc'st quarter, precipitation of 
driest quarter, precipitation of warmest quarter, temperature 
annual range (Fig. SI; [51], http://worldclim.org). The contem- 
porary (~ 1,950-2,000 AD) bioclimatic variables were selected 
based on their ability to capture seasonal trends and species- 
specific differences between M. gentilis and M. onohychidis. 
Differences in the distribution of the study species for each 
bioclimatic variable were tested with the WUcoxon rank sum test 
with significance (P) based on an alpha of 0.95 and visualized with 
violin plots (Fig. SI). Data concentration ellipses were estimated 
using the k-means clustering method (K= 2) to visualize contem- 
porary record assignment, as well as assign the environmental 
space of the La Brea Tar Pits during the LGM to either M. gentilis 
or M. onobiychidis. 

Geographic distribution of habitat suitability. To esti- 
mate the geographic distribution of habitat suitabilit)' in the 
western U.S. at present and during the LGM (-21,000 y BP), 
habitat suitability models (HSM) were constructed under the 
principle of maximum entropy with MaxEnt v3.3.1 [52]. Each 
species' HSM was constructc-d from the specimen locality records 
and bioclimatic variables selected for the principal component and 
cluster analysis (see Environmental Space Use). Under the 
principle of niche conservatism [53], we estimated the geographic 
distribution of habitat suitability of both species during the LGM 
with analog bioclimatic variables of their contemporary biocli- 
matic niche. To provide a robust estimate of LGM habitat 
suitability, we allowed for the HSM to exceed contemporary 
bioclimatic extremes in the event that maximum entropy was 
achieved during model construction. Constructing LGM HSMs 
that are younger than the fossil records of the study species 
guarantees that both M. gentilu and M. onobrychidis, or a common 
Litomegachile ancestor were already in existence. The spatial 
resolution of all bioclimatic variables used in the models is 2.5 
arc-minutes. Data processing and geographic visualization of the 
HSM were implemented in ArcGIS 10.1 (ESRI, Redlands, CA). 
Model performance was evaluated in MaxEnt using the area 
under curve (AUC) statistic: and a jackknife of variable importance 
(Table SI). The AUC is constrained between zero and one where 
values closer to one suggest better model performance. Models 
were averaged over 50 repUcates, with 10-fold cross validation for 
each replicate for model validation. 

Statistical analysis. Outside of MaxEnt, all analyses were 
conducted in R v3.3.0 [54]. Principal components were visualized 
with the scatterplot function in the car library. The k-means 



clustering method was performed with the kmeans function in the 
cluster library. Violin plots were visualized with the vioplot library. 

Supporting Information 

Figure SI Violin plots of eight bioclimatic variables 
associated with the distributions of M. gentilis and M. 
onobrychidis. The width along each \iolin plot represents the 
frequency of specimen records associated with a measurement of 
the bioclimatic variable under observation. Wider widths reflect a 
higher frequency of specimen records, whereas thinner widths 
reflect a lower frequency of specimen records. Two-sample 
Wilcoxon tests were performed to test for precipitation/temper- 
ature differences between the species for each bioclimatic variable. 
Alpha levels for each test were set at alpha = 0.95. /'<0.05 suggests 
significant differences in bioclimatic space use between the two 
species. M. gen = M. gentilis and M. ono = M. onobrychidis. 
(TIFF) 

Figure S2 Contemporary habitat suitability distribu- 
tion of M. gentilis and M. onobrychidis across an 
elevation gradient in southern California. Kendall T rank 
correlation coefficient estimates for M. gentilis revealed a positive 
and significant correlation between HS and Elevation (T = 0.31, 
P<0.05), and a negative and non-significant correlation for M. 
onobrychidis (T = 0.31, P=0.52). 
(TIFF) 

Table SI Analysis of bioclimatic variable importance 
for the MaxEnt habitat suitability models (HSM) of M. 

gentilis and M. onobrychidis. Values represent averages over 

50 replicate runs of each species' habitat suitability model (HSM). 

(PDF) 

Dataset SI Natural history records of M. gentilis and 
M. onobrychidis utilized to construct environmental 
space use and habitat suitability models (HSM). AU 

specimens of ea[:h species present in the U.S. National Pollinating 
Insect Collection Database are presented in the tabk;. Megachile 
hrevis onobrychidis has been synonymized as M. onohiyihidis. The 
former name is retained in the data table. Duplicate specimen 
records for each species per unique locality were removed for the 
final analyses. 
(XLSX) 

Video SI LACMRLP 388Ea. Rotating video nest of 
opaque to transparent nest cell containing male pupa. 

(MP4) 

Video S2 LACMRLP 388Eb. Rotating video nest of 
opaque to transparent nest cell containing female pupa. 

(MP4) 

Video S3 LACMRLP 388Ea. Rotating dorsal to ventral 
video of male pupa. 

(MP4) 

Video S4 LACMRLP 388Ea. Rotating dorsal to ventral 
video of female pupa. 

(MP4) 

Video S5 LACMRLP 388Ea. Rotating lateral video of 
male pupa. 

(MP4) 

Video S6 LACMRLP 388Eb Rotating lateral video of 
female pupa. 

(MP4) 
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Video S7 LACMRLP 388£a. Rotating quarter edge of 
nest cell. This video shows cross-sections of leaves from different 
angles. 
(MP4) 
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